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ABSTRACT 

Gravitational lensing is potentially able to observe mass-selected halos, and to measure 
the projected cluster mass function. An optimal mass-selection requires a quantitative 
understanding of the noise bahavior in mass maps. This paper is an analysis of the 
noise properties in mass maps reconstructed using a maximum likelihood method. 

The first part of this work is the derivation of the noise power spectrum and 
the mass error bars as a straightforward extension of the Kaiser & Squires (1993) 
algorithm to the case of a correlated noise. A very good agreement is found between 
these calculations and the noise properties observed in the mass reconstructions limited 
to non-critical cluster of galaxies. It demonstrates that Kaiser & Squires (1993) and 
maximum likelihood methods have similar noise properties and that the weak lensing 
approximation is valid for describing it. 

In a second stage I show that the statistic of peaks in the noise follows accurately 
the peak statistics of a two-dimensional Gaussian random field (using the BBKS tech- 
nics) if the smoothing aperture contains enough galaxies. This analysis provides a 
full procedure to derive the significance of any convergence peak as a function of its 
amplitude and its profile. 

It is demonstrated that, to a very good approximation, a mass map is the sum 
of the lensing signal plus a 2D gaussian random noise, which means that a detailled 
quantitative analysis of the structures in mass maps can be done. A straightforward 
application is the measurement of the projected mass function in wide field lensing 
surveys, down to small mass over-densities which are individually undetectable. 
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• i-H 1 INTRODUCTION ticular a joint analysis of clusters of galaxies from weak lens- 

J\. r 1 1 ingi Sunyaev-Zeldovich, X-rays and galaxy catalogues would 

, • Non-parametnc mass reconstruction from gravitational shear , i 'i i a i . i • ■ . i -i i ■ , r i i 

M , , , , be valuable. A related issue is the possible existence oi dark 

cn ' IS now a very well understood problem, and was widely re- , /tt , , ■ i r.r.r,\ i i i , , i , , ■ r 

. . . 1 „ , , , ^ ^ . ' , . , . tenses (Hattori 1997), and already, the detection or a mass 



clump not obviously associated with light (Erben et al. 1999) 
raised the question of quantifying accurately the significance 



fined over the last few years. It is no longer a technical is- 
sue to reconstruct clusters (see Squires & Kaiser 1996 for 

a detailed review) as well as large scale structures (Van „, , , 

OI il3/iO ClStSCtlOIlS 

Waerbeke et al. 1999). Although the computational speed „, i . i. r,^ ■ ■ 

' Ihe understanding or the noise properties m mass maps 

or mass reconstruction might be too low lor large scale sur- ■ r r i . . . i 

u J- c ID • innnN i,- • • theretore necessary tor anyone who wants to measure the 



veys (Lombardi & Bertin 1999), this is not a major concern 
yet according to the actual size of the lensing surveys. 

Naturally a next step in gravitational lensing is the 



mass distribution from gravitational lensing in order to con- 
strain cosmological models. This point was understood by 

Schneider (1996) who constructed the so-called S-statistic, 
search lor a better use ot the information present m the ■ i ^ ^-c ■ -n r i tt • i ^ 

aimed to quantity the significance ot mass peaks. Using boot- 
mass maps. Assuming that all the problems concerning the , . r.ii.ii - r r i 

strapping ot the data, he derives the frequency ot observed 

residual of Point Spread Function corrections are or will be 



mass peaks. The S-statistic is similar to the the signal-to- 
noise derived in Kaiser & Squires (1993) (and confirmed 
later in numerical simulations by Wilson et al. 1996) but it 

has the interesting feature that works directly with the ellip- 
(Reblmsky & Bartelmann, 1999) as a probe ol cosmology i- r i • ^ ^ ^- ■ ■ ^\ 

,„ , , „ , , „ , , , , N , ticity ot the galaxies (no mass reconstruction is required), 

(Oukbir & Blanchard 1997, Rebhnsky et al. 1999). In par- J & v h 



solved, there is useful information hidden in the noise of 
mass maps that can be extracted. An important application 
of lensing is the possibility to built mass-selected samples 
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and the smoothing filter can be optimized in order to in- 
crease the signal-to-noise. 

So far, the noise analysis and the statistical significance 
of mass reconstructions have not been studied as much as 
the mass reconstruction methods themselves. We have only 
some global clues of the behavior of the noise, which above 
all, seem to depend on the mass reconstruction method 
(Squires & Kaiser 1996). Seitz & Schneider (1996) identi- 
fied several characteristics of the noise in their numerical 
simulations: due to the smoothing, the power spectrum is 
not flat and driven by the intrinsic ellipticity of the galax- 
ies, their number density and the smoothing window scale. 
Lombardi & Berlin (1998a, 1998b) gave the first analyti- 
cal estimate of the variance of the total mass of a recon- 
structed cluster, taking into account the noise correlation 
induced by the smoothing window. They have found differ- 
ences with the numerical results of Seitz & Schneider (1996) 
that they attributed to technical differences. They worked 
with shear line-integral reconstruction methods, and found 
that the noise is minimized when curl-free kernels are used. 
As discussed in Seitz & Schneider (1996) the curl- free ker- 
nels have a noise amplitude similar to the Kaiser & Squires 
(1993) (hereafter KS) reconstruction, but the line- integral 
reconstruction is a non-local process, which complicates the 
noise analysis (Lombardi & Berlin (1998a, 1998b)). On the 
other hand, the Maximum Likelihood (hereafter ML) algo- 
rithm (Bartelmann et al. 1996) performs a local reconstruc- 
tion, but its noise properties can hardly be studied analyt- 
ically. In fact, the KS scheme seems to be the easiest ap- 
proach to an analytical description of the noise. 

It has been demonstrated in the case of lensing by large 
scale structures (see the simulations performed in Van Waer- 
beke et al. 1999) that the noise in ML reconstructions is ac- 
curately described by the weak lensing approximation (and 
therefore by the KS equations). Continuing this work, this 
paper is a study of the noise properties of ML methods in 
the case of non-critical clusters of galaxies with correlated 
noise, when the smoothing is taken into account. 

This paper shows that the noise in ML reconstruction 
agrees very well with the analytical prediction made from 
the weak lensing approximation. It demonstrates that the 
mass variance derived in Lombardi & Berlin (1998b) is valid, 
and, to a larger extent, that the noise is uncorrelated with 
the signal even for non-critical clusters of galaxies. This is 
also consistent with the empirical analysis done in Seitz & 
Schneider (1996). Assuming that the smoothed intrinsic el- 
lipticities are described by the statistics of a two-dimensional 
Gaussian random field, the BBKS technique (Bardeen et al. 
1986) is used in order to quantify the significance of the over- 
densities in ML mass reconstructions. This analysis provides 
a quantitative way to identify the structures in mass maps 
and to extract information from the noise: this result is ver- 
ified on numerical simulations, which shows that to a very 
good approximation a mass map is the sum of a completely 
known Gaussian field plus the lensing signal. 

Section 2 gives some general considerations about lens- 
ing theory used in this work. The noise model is developed 
in Section 3, where each noise contribution is identified. An 
estimator for the variance of the mass in a given area is 
calculated. In Section 4, the noise model and numerical sim- 
ulations of cluster lensing are compared. In Section 5 the 
significance level of a reconstructed mass halo based on a 



2D Gaussian field statistic is derived, which is then checked 
against numerical simulations. 



2 GENERALITIES 

We consider a lens of surface mass density E(0) located at a 
redshift zi along the line of sight 9, and a single source plane 
is put at a redshift Zg. At first approximation, the Universe 
is equivalent to a single lens, so there is no loss of generality 
here to consider a single lens. The projected gravitational 
potential <E> and the surface mass density of the lens are 
given by the 2-dimensional Poisson equation 

A$(0) = 47rGE(0). (1) 

The critical density (for which the lens produces multiple 



images) is given by 
c' 0„., 



(2) 



where T>ij is the angular diameter distance between the red- 
shifts Zi and Zj. For convenience, the lensing potential is 
defined as </> = ^ / {2nGT,crit)- The ratio 



«(0) = 



E(0) 



(3) 



is the convergence of the lens, which describes the isotropic 
deformation of a light beam coming from the sources. The 
light beam is also sheared, which modifies the shape of the 
source galaxies. The shape of a galaxy is described by the 
second moments of its centered surface brightness profile 
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(4) 



The source and image shape matrices are related by the 
equation 



(5) 



where A is the magnification matrix which depends on the 
convergence k and the shear 'y of the lens. 
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K and "f are fully determined by the second derivatives of 
the lensing potential via 



K = A(t>/2- 71 
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The ellipticity is defined in the complex plane as 



1 - 2itf 
1 + r 



(7) 
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where r is the axis ratio of a galaxy and i9 the orientation. 
Eq.(|^ written in terms of source and image ellipticities is 



l-hg*e 



(9) 



where g = 7/(1 — k) is the complex reduced shear. The 
observable is <5 = 2g/(l + jgp), but this paper is restricted 
to non-critical lenses, for which \8\ < 1, so that g is an 
observable as well. In that case g can be immediately read 
off from the galaxy ellipticities; as shown by Schramm & 



Kayser (1995) and Schneider & Seitz (1995) the mean ob- 
served ellipticities of the galaxies is an unbiased estimate of 
the reduced shear: 



(10) 



The mass reconstructed from shear measurement alone is 
not fully constrained because any transformation of the lens- 
ing potential of the form 



(11) 



leaves the reduced shear g unchanged. This is the so-called 
mass-sheet degeneracy: the mass is determined up to an ad- 
ditive constant and an appropriate scaling (Eq.(pT|) implies 

K — > 1 — A + 



3 NOISE CORRELATION AND MASS ERROR 
BAR 

3.1 Weak lensing approximation 

As already mentioned, the weak lensing approximation lin- 
earises the lensing equation and dramatically simplifies the 
noise analysis. The weak lensing limit of Eq.(b|) is 



JS) 



7 + (12) 

Thus (e'^') = 7 becomes an unbiased estimate of the shear. 
In the Fourier space Eq.(|7|) takes the form (Kaiser 1998): 



K(k) = Ca{k)jcO^), 



(13) 



where there is an implicit summation over a = (1,2). 
Cq — {cos{2(fi)) , sm{2ip)) , and k = k{cosLp,smip). An un- 
biased (but extremely noisy) estimate of the convergence n 
is Kn(d) = i^(d) + n{6), where n{6) is the noise part, which 
depends on the intrinsic ellipticities of the galaxies and their 
spatial distribution. The Fourier transform of k„ is 



ft„(k) = c^e'^JHk) = k{k)+c^ei^\k). 



(14) 



In order to avoid the infinite variance problem (caused by the 
discreteness distribution of the galaxies), and to have a mass 
map with a reasonable amplitude of noise, the convergence 
is reconstructed from the smoothed galaxy ellipticities, so 
smoothed fields are used in the following. Uppercase letters 
refer to smoothed fields while lowercase letters refer to un- 
smoothed fields (e.g. K for the smoothed k, T for "f and S 
for e). The smoothing function is denoted by Vl^(|0j), the 
total number of galaxies in the field by Ng, and the total 
field area by 5. Galaxies are at positions 6i. The recon- 
structed mass map is now the smoothed convergence with 
an unbiased estimate given by 
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e,\)n{e. 



(15) 



The number density of galaxies is Ug — Ng/S. It is worth 
noting that Eq.(p^) implies that the reconstructed conver- 
gence of a smoothed shear field is identical to the initial 
smoothed convergence. In other words, smoothing and re- 
construction are commutative processes in the frame of the 
weak lensing approximation. 



3.2 Noise correlation function 

The reconstructed convergence Km(0) is the sum of the true 
convergence K(6) and a correlated noise N{e). If S<^' and 
S*-^' are the smoothed e*^' and e''^' then 

J<iv(6») = K{e)+N{e) = y" dk e^*-^c„(fc)Ei^'(k). (16) 

If the smoothing length is small enough (arcmin scale) the 
signal is approximately constant within the aperture. This 
is certainly not true for cluster mass profiles, but it should 
still give a reasonable description of the noise. It implies 
that the sum over the galaxy ellipticities within the aperture 
is equal to the smoothed shear plus the intrinsic ellipticity 
contribution (which is Eq.(|l3) for the smoothed fields): 



T,^'\e) ~ r(0) + — y VK(0 - e,)e^^\ei 



(17) 



Inserting this into Eq.(^^ gives the analytic expression of 
the noise term: 



KN{e)-K{e)~^Y^ / dkiy(k)e"* <^'-^' 



Cc(fc)e<.(6»0,(18) 



where W^(k) is the Fourier transform of the window func- 
tion. The ensemble average of the square of Eq.(p^) is the 
noise variance (N'^{9)). It is calculated by averaging over the 
intrinsic ellipticities and the galaxy positions. Let A be the 
operator which averages over the intrinsic ellipticities and P 
the operator which averages over the galaxy positions. The 
former is defined as 



A 



-6'). 



(19) 



is the Kronecker symbol, and CTe the intrinsic ellipticity 
dispersion. This is a consequence of the fact that the galaxy 
ellipticities are assumed to be spatially uncorrelated. The 
galaxy position average operator P depends on the angular 
clustering of the galaxies within the aperture (Peebles 1980): 

P. / ^^4^i^i + ^^\e,-e,\) + ...) (20) 



where S is total field area, and a;(|02 — 0i|) is the angular two 
point correlation function of the source galaxies. This clus- 
tering term is generally ignored, and accounts for possible 
increase of the variance of N{0) in the case of non-uniform 
sampling of the source galaxies within the aperture. Its con- 
tribution can be potentially important if the source galaxies 
belongs to a thin redshift slice (as the angular clustering 
approaches the spatial clustering). Connolly et al. (1998) 
measured on the Hubble Deep Field the angular correlation 
function of galaxies in the redshift interval [1, 1.2], which cor- 
responds to the most likely source redshift for weak lensing 
surveys. They found an angular correlation length of 3 arc- 
sec, and therefore negligible for our purpose. In practice the 
redshift slice of the lensed galaxies is much larger [0.5, 1.2], 
which makes the contribution of the angular clustering to- 
tally negligible (Le Fevre et al. 1996) . As a consequence, the 
ensemble average of the local variance of the convergence is 
simply given by: 

,2 



V[K[N\e)\] 
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dk|iy(k)|^ 



(21) 



as derived in Kaiser & Squires 1993. The noise correlation 
function is then 



p[A [N{e)N{e')]] 



2n„ 



dk e*'^'-^)|VK(k)l^ 



(22) 




3.3 Average mass error bar 

The average mass error bar AM in an aperture (the to- 
tal mass of a cluster for instance, or the mass of a sub- 
structure) is now calculated. Provided that the redshifts 
of the background galaxies are all the same, the projected 
mass M(6) in the direction 6 is M{9) — T^critizi, Zs)K{6), 
where T,crit{zi, Zs) is the critical surface mass density of the 
lens which depends on the lens configuration (Eq.(|^). If 
the background galaxies were distributed in redshift, the re- 
duced shear becomes redshift dependent but it can be easily 
accounted for (Seitz & Schneider 1997), and with no loss of 
generality the analysis is restricted to a single source red- 
shift. 

The variance of the total projected mass estimated 
within an area 5 is by definition 

AM^ = f f dO dO'P [A [Nie)N{e')]] . (23) 

Js Js 

According to Eq.([22[) the variance is finally 
-2 r 



2 n-g j_j 



dr / dk e* "'|M/(k) 



(24) 



3.4 Gaussian smoothing 

The Gaussian window has a width cr, and is normalized: 



(25) 



Inserting this into Eq.(^2) gives the noise correlation func- 
tion 



{N{e)N{e')) = 
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(26) 



The variance defined by Eq.(|2^) for a square area S = 
[2a, 2a] is then straightforward to calculate; 



AM^ 



4a* 



exp 



^erf (^^\ 
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(27) 



3.5 Top-hat smoothing 

The normalised top-hat of radius a is 

W^d^i) = \ if|01 <a 

= otherwise 
The noise correlation function is then 



(28) 




Figure 1. The cluster mass model used in this work sampled on a 
42x42 pixels grid, with a maximum convergence Kmax = 0.3. The 
lens, defined by Eg. (^l]) has a redshift z = 0.2, and an Einstein-de- 
Sitter Universe is assumed. This is the mass distribution caused 
by two truncated isothermal potential with a core radius of 5 
pixels. 
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= otherwise 



if|0' -e\<2a 
(29) 



where s — 0|) is the overlapping area of two identical 
circles of radius a distant hy A — \0' — 6\ 
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4 CLUSTER LENS SIMULATIONS 
4.1 Description 

The description of the noise given above is now compared 
to mass reconstruction simulations using the lens compara- 
ble in lensing strength to a non-critical cluster of galaxies 
{i^max ~ 0.3). The lens potential <j!> is the sum of two non- 
singular isothermal spheres with a core radius Rc, located 
at Oi and 02: 



C2 Vos 
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\e-92 
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(31) 



The velocity dispersion is a^c ~ 900 km/s, Rc = 5, 
and we assume a lens redshift zi — 0.2. The distances are 
calculated in an Einstein-de-Sitter Universe. Figure |^ shows 
the convergence k = A'i/) of the lens. The shear caused by 
the lens is applied to a background galaxy population at 
redshift Zs = 1, with an ellipticity distribution given by 



p{e^^^) cx exp 
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(32) 



with (Te = 0.12 (corresponding to a typical axis ratio of 
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Figure 2. The noise correlation function measured in the sim- 
ulations (thin lines), compared to the model (thick lines). The 
dashed lines correspond to the Gaussian smoothing and the solid 
lines for top-hat smoothing. Error bars are calculated from 50 dif- 
ferent simulation realizations, they are displayed for the gaussian 
smoothing only, but they are similar for the top-hat case. 

0.8). 1000 background galaxies are randomly placed on a 
42 X 42 pixels grid. An observed simulated shear map is then 
calculated, and the lens mass distribution is reconstructed 
with a maximum-likelihood algorithm (Bartelmann et al. 
1996). The maximum likelihood reconstruction minimizes 
the function defined as 

^'^\s{'t> guess) - Sobsl'^ , (33) 
id 

where (pguess is a guessed lensing potential. The reduced 
shear g{(j>guess) is calculated from (pguess, and gobs is the ob- 
served reduced shear. The reduced shear is calculated from 
the finite difference scheme described in Van Waerbeke et al. 
(1999). The comparisons are done between the reconstructed 
mass map and the initial mass smoothed with the same win- 
dow used to smooth the shear for the reconstructions (mass 
reconstruction and smoothing should commute according to 
the noise model) . For each type of smoothing window (top- 
hat and Gaussian) 50 simulations are performed. The widths 
of the smoothing windows are a = 3 pixels, therefore there 
is a mean number of 16 galaxies per aperture. An edge of 
width 3 pixels is systematically removed from the analysis 
because of its intrinsic higher noise due to the change of the 
finite difference scheme at the boundaries of the grid and 
the smaller number of galaxies in the aperture. The cali- 
bration constant A (see Section 4.2) is calculated by forcing 
the total mass of each reconstructed map to be equal to the 
total mass of the smoothed mass distribution of the lens. 
Therefore the variance of the total reconstructed mass map 
is zero, so the analysis is restricted to half of the grid size. 
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Figure 3. The crosses on the first plot show the total mass 
Mrecon measured in the simulation for different sizes of the aper- 
ture S and for the gaussian smoothing case. For the smallest 
aperture, which corresponds to the pixel size, Mrecon — 0.1. The 
largest aperture has a width of 18 pixels (about half of the sim- 
ulation grid), and Mrecon — 36. The vertical axis displays the 
smoothed mass distribution of the model. The solid line shows 
a perfect correlation between Mrecon and Mmodel- The recon- 
structed mass is overestimated by less than 2%. The diamonds 
on the second plot show the standard deviation of the mass for 
different aperture sizes. The horizontal axis corresponds to the 
measure of AAf in the simulation, and the vertical axis is the 
prediction from Eg. (p?]). The solid line shows a perfect correla- 
tion between AMrecon and AMmodei- 

4.2 Effect of the mass-sheet degeneracy 

For each reconstruction, the mass map is rescaled in such a 
way that the total mass measured on each realization is the 
same as the total mass in the smoothed lens model. 

Such a rescaling changes the amplitude of the noise cor- 
relation function: if Kreai is the real mass distribution, the 
reconstructed mass map Krccon is related to the real mass 
map by Krccon = 1 — A -|- \Kreai + N , where A is the mass- 
sheet degeneracy constant and A*' the noise term. Therefore 
the estimate of the real mass distribution is Krcai + -r- where 
the noise scales now with A .It follows that the noise cor- 
relation function of the reconstructed mass map scales as 
A~^, according to Eq.(p2[). In that sense, error bars of K de- 
pend on the absolute value of K, although it is not explicit 
in Eq.diil), and this rescaling has to be taken into account 
when the amplitude of the noise correlation function is cal- 
culated. 



4.3 Correlation functions 

Figure ^ is a comparison of the measured noise correla- 
tion function to the predictions (Eq.(|26|) for the Gaussian 
smoothing and Eq.(G9|) for the top- hat smoothing). The 
dashed lines correspond to the Gaussian smoothing and the 
solid lines to the top-hat. The thin lines display the measure- 
ment in the simulations, they are in very good agreement 
with the predictions (thick lines). As expected, the correla- 
tion drops to zero for the top-hat window when |r| = 2a. 
The amplitude of the correlation function is twice the value 
for the top-hat compared to the Gaussian, as expected from 
Eq.(|26|) and (|i). 



4.4 Mass error bar 

The average variance of the mass in an area S = 2a x 2a 
is given by Eg. (p^) . This quantity does not depend on the 
signal according to our hypothesis. As shown by the plot 
at the top of Figure ^ the reconstructed mass is an unbi- 
ased estimate of the smoothed initial mass distribution (as 
expected from Eq.(|l^)). The reconstructed mass is slightly 
overestimated by less than 2%. 

The diamonds on the bottom plot of Figure |^ shows the 
error Eq.(|^ {AMmodei), compared to the measured error 
on the simulations (AMrecon)- The different diamonds cor- 
respond to different sizes of the smoothing area S = 2a x 2a 
{Arecon — 0.1 for 2a — 1 pixel, and Arecon — 1.5 for 
2a — 18 pixels). The agreement between the model and 
the measurement in the simulation is good, however, there 
is a tendency of the measured variance to be slightly un- 
derestimated for 2a ~ 18. This is a consequence of the nor- 
malisation of the mass maps (the total measured mass is 
explicitely set to to the total mass in the model, so the mea- 
sured variance is zero at the scale 2a = 42 pixels). The errors 
are independent of the position of 5 on the grid, which con- 
firms that signal and noise are uncorrelated to a very good 
approximation. 




Figure 4. One realization of mass reconstruction with no lensing 
inside. All the structures here are noise. The high density peak 
at the center is a 3cr fluctuation (k = 0.081). The dark square 
represents the subfield where the statistical analysis is done (in 
order to avoid the noisier boundaries), and the circle at the top 
right has a 0.66' diameter, representing the size of the Gaussian 
window width. The field is 14' X 14' large, with a number density 
of 30 galaxies per arcmin,'^. The intrinsic ellipticity dispersion is 
(Tj = 0.17. 



I propose here to use the statistics of 2D Gaussian ran- 
dom fields as a tool to analyze the weak lensing mass maps. 
With this formalism, the significance and the morphology of 
the structures can be evaluated analytically and compared 
wih the data, which is a way to go beyond the detection 
of dark halos done in Schneider (1996), based on a local 
signal-to-noise analysis. 



5 SIGNIFICANCE OF PEAKS 

The reconstructed mass maps are generally very noisy, even 
with an arcmin scale smoothing, and it is difficult to dis- 
entangle real and noise over-densities. Figure ^ is an exam- 
ple of mass reconstruction when no lenses are present. It 
shows fiuctuations which can be misinterpreted as real mass 
over-densities. It makes particularily difficult the selection 
of objects based on their mass, and in practice we use a 
high signal-to-noise cut-off to be sure to minimize the num- 
ber of spurious peak detections (Schneider 1996, Reblinsky 

6 Bartelmann 1999). According to Section 3.2, the noise in 
the reconstructed mass maps is equivalent to the convolution 
of a point process by a smoothing aperture, therefore it is 
expected to approach a Gaussian field distribution when the 
aperture is larger than the mean angular distance between 
galaxies. The nice feature of Gaussian fields is that they are 
fairly well understood (Bardeen et al. 1986), in particular 
the statistics of two-dimensional Gaussian random fields as 
been worked out in detail for the case of cosmic microwave 
background fluctuations (Bond & Efstathiou 1987, hereafter 
BE). 



5.1 Height of peaks 

Let N{6), the noise part of Eq.(|l|), be a 2D Gaussian ran- 
dom field. The field and its derivatives (V = [N,N',N"]) 
are distributed as 
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where the correlation matrix M is 
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The cr„ are determined from the noise power spectrum 

^k'"mk)\'), (36) 

which can be calculated from the two-point correlation func- 
tion: 
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The zero-lag value of the correlation function is ao = 
(iV(O)iV(O)). The marginalisation of Eq.(|3) over the field 
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Figure 5. In all the plots, the dashed curve is the probability distribution function that a randomly chosen pixel has a value uag 
(Eq.disI)), and the dot-dashed line is the expected probability distribution of peaks extrema, according to Eg. (|4o[). The solid line in the 
top-left plot is the measured distribution of all the pixels in the 90 simulations (described in Section 5.1.) limited to the sub-field defined 
by the dark square as explained in Figure ^ The solid line in the bottom-left plot shows the distribution of extremas measured in the 
same sub-fields of the simulations. In the top-right plot, the solid line is the measured distribution of the value of all the pixels in the 
whole field, and in the bottom-right plot, the solid line is distribution of extrema measured in the whole field. 



derivatives gives the -Gaussian- probability that any pixel 
randomly chosen in the field has the value vao: 
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In order to describe the probabihty distribution function of 
peaks, the following parameters are introduced: 7 and 9^ 
are respectively the width of the power spectrum and the 
characteristic angular scale of the field: 



7 = 
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A peak occurs, by definition, where the gradient of the field 
vanishes {VgN = 0), and it has a given height i/ctq. Its prob- 
ability distribution function 'P(//|peafc) can be calculated 
from the number density of minima and maxima Nmax{i') 
and Nmin (v) (see BE for the details) : 



'P{v\peak) — 



where 
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Nmax{v) ^N^^n{-u) = exp {~^'' l^) ^^J^ . (41) 



The function G{'y, x) depends on the variable = 2(1 — 7^ 
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The denominator of Eq. ( [40| ) has an analytical solution which 
can be derived straightforwardly with a formal calculator, 
but it is not worth to give its expression here. 

For a Gaussian smoothing (Eq.(^5|)) the spectral pa- 
rameters take the values 7 = l/\/2 and 9^ = a/\/2, and 
therefore Vivlpeak) only depends on p. The denominator of 
Eq.(Eol) is: 
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(43) 



It means that once the mass map is normalised by ctq the 
probability distribution of peaks is always the same, inde- 
pendently of the smoothing scale, the number density of 
background galaxies, and the ellipticity dispersion a^. 

The analytical predictions of Eq.(|38[) and Eq.(^) are 
now compared to numerical simulations. Ninety mass recon- 
structions have been performed, with no lensing signal in- 



side. The total field of view for each realization is 4100 x 4100 
pixels (equivalent to a 14' x 14' area at the Canada France 
Hawaii Telescope), and it contains about 30 gal/arcmin . 
The intrinsic ellipticity dispersion is = 0.17 (using the 
distribution Eq.(H)), and the smoothing aperture is given 
by Eq.(p5|) with a = 0.33', which corresponds to an average 
of 10 galaxies per aperture. The "mass" field is reconstructed 
using the maximum-likelihood algorithm. Since the aperture 
is truncated at the edge of the field, only a sub-set of pix- 
els in the simulations are used, those enclosed in the square 
area defined on Figure ^. This Figure also shows the scale 
of the smoothing window compared to the whole field. 

Figure |5| shows the measured probability distribution 
functions in the field and for peaks, for the whole field of the 
simulations, and for the sub-field as described in Figure ^ 
The measurements are represented by the solid lines, while 
the dashed and dot-dashed curves show respectively the un- 
constrained (Eq.(^8|)) and the constrained (Eq.(pc|)) distri- 
butions. For the centerred sub-field it is remarkable that the 
peak distribution prediction fits so well the measurements. 
It shows that the non-linear ML reconstruction give a noise 
comparable to the noise in weak lensing approximation, con- 
firming our previous results. However it is clear that when 
the measurement is done in the whole field, the statistic is 
not gaussian. This is an effect caused by the boundary of 
the reconstruction: at the edge of the field, the smoothing 
aperture is truncated and it contains less galaxies than in 
the central part of the simulations. Therefore the peak anal- 
ysis has to be restricted to the area where the aperture is 
not too close to the edge. Two other sets of 20 simulations 
each, one with — 0.4 and one with a uniform distribution 
of ellipticities between -0.2 and 0.2 (instead of Eq.(^) give 
the same result, which shows that the Gaussian behavior of 
the noise does not depend on the details of the ellipticity 
distribution. The distribution of peaks is not Gaussian, as 
expected for Gaussian fields. It shows that simple signal- 
to-noise analysis based on the local variance with Gaussian 
probabilities is wrong. For instance a value at 3a in the field, 
with Viiy) — 0.004 corresponds to Vivlpeak) = 0.03 for a 
peak, almost a factor of ten higher. 

Because of the discreteness of the spatial distribution 
of the galaxies a Poisson-like noise is expected rather than 
a Gaussian noise. The fact that the noise is not Poisson is 
a consequence of the central limit theorem which makes the 
distribution more Gaussian because there is an average of 
10 galaxies per aperture, which is much larger than 1. 

The peak statistic described here is very easy to im- 
plement and it requires only the calculation of cto. In other 
words, only the galaxy number density, the variance of their 
ellipticities and the shape of the smoothing window are 
needed. 



5.2 Profile of peaks 

It is useful to add the morphological information to the peak 
analysis in order to constrain the mass distribution better. 
Let Vishapelv) be the probability that a peak of height i/ao 
has a given shape, and 'P{shape,h') the global peak proba- 
bility written as: 



'P{shape,v) = 'P{iy\peak)'P{shape\i') 



(44) 



BE give the analytical formulae describing the mean profile 
of a peak and its dispersion. This calculation depends on the 
averaged curvature of a peak of height i/. The curvature x is 
defined from the eigenvalues (Ai, A2) of the second derivative 
matrix of the noise dijN: 

X- = (Ai + A2)/o-2. (45) 

Following BE, all distances and derivatives are calculated 
with respect to the reduced coordinates UJ = r/di,. The re- 
duced correlation function is defined as tjj{uj) = {N{9)N{9+ 
Lj)} / {N{0)N{0)) . The profile around a peak, averaged over 
the curvature and the ellipticity, is Gaussian distributed 
with a mean {N{u:)\u) and a dispersion {N'^{u))\u) (Eqs. 
A2.9 and A2.8 of BE): 
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For 2D Gaussian fields, the conditional mean value of the 
average curvature {x\v) can be calculated analytically: 
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(47) 



G(7,7!/) ■ 

The function H{'y,x)/G{'y,x) tends to zero for high peaks 
(u — > 00). G{'y,x) was given in Eq.(|42[) and H{'y,x) is: 



(48) 



The quantity quantifies how any observed peak profile of 
height vao Nobs{d, v) is different from its expectation value: 



(49) 



This follows a statistic with the number of degrees of 
freedom n given by the number of positions u: involved in the 
summation. Thus the probability 'P{shape\v) that a peak of 
height 1/(70 having a definite shape is a noise fluctuation is. 



(50) 



The probability Vlshape^v) that a peak of height vao with 
a given shape is a noise fluctuation is then given by Eq.(w4), 
which contains the probabilities ( |4o| ) and (|5o|). Figure ^is 
an example of mean profile and its dispersion derived from 
Eqs. 



6 CONCLUSION 

It is possible to understand the noise properties of recon- 
structed mass maps with rather simple arguments. The noise 
in ML mass reconstruction is very close to a 2D Gaussian 



A 




Figure 6. Mean profile around a peak of height Tctq (top 
curve) and l.Scro (bottom curve) with the noise correlation 
function given by Eq. (2^. The parameters are ct = 3, rig = 
0.6 galaxies/pixel^ , and = 0.12. The dashed curves represent 
the it 1(7 deviations. 



random field, which iias several implications for the analysis 
of mass maps: for cluster lensing, the "blurred" area around 
cluster of galaxies where the identification of mass clumps 
and substructures is uncertain can in fact be described ac- 
curately. It is possible to give the probability that the fluc- 
tuations are real mass-over-densities, and to give their mass 
with an error bar. For large scale surveys, small mass concen- 
trations (not visible individually) can be detected statisti- 
cally as a deviation from the expected noise Gaussian statis- 
tic. Groups and small clusters should show up in the statis- 
tical distribution of peaks. Preliminary results using ray- 
tracing simulations (Jain, Seljak, White 1999, Jain private 
communication) already confirmed this, but the detailled 
analysis of peak statistic in realistic cosmological models is 
left for another work. 

The noise autocorrelation function and the mass error 
bars were derived and compared to numerical simulations of 
non-critical cluster lenses. The agreement between the pre- 
dictions and the simulations proves that the noise, caused 
by the intrinsic ellipticities and the discrete distribution of 
the galaxies, is nearly uncorrelated with the lensing signal 
despite the non-linear lensing equation. Consequently, the 
weak lensing approximation applies, which simplifies con- 
siderably the noise analysis. It is worth to remind that this 
is true for non-critical clusters only. However it is reasonable 
to think that the smearing of the signal due to the smoothing 
implies that the approximation still works for noise analysis 
in critical clusters (if the smoothed mass distribution be- 
comes non-critical). 

A BBKS-type approach was used to predict the statis- 
tics of noise peaks in mass maps. Again the numerical sim- 
ulations are accurately fitted by these predictions. The fact 
that the noise approaches a Gaussian distribution is due to 
the large mean number of galaxies in each aperture (about 
10). There is certainly a limitation of the method used here: 
if the smoothing window is too small or too "fuzzy" it should 
deviate from a Gaussian statistic. This departure from Gaus- 



sianity was observed in the S-statistic of Schneider (1996) 
because he used a compensated filter. A work worth to be 
done is to investigate the limits of the method, and even- 
tually to consider the Poisson behavior of counts in small 
cells, as in Szapudi & Colombi (1996). However, in the con- 
text discussed here, it is reliable to use Gaussian statistics. 

I discussed so far only the ML mass reconstruction 
method, and compared it to the analytical KS. Can this 
analysis be used for other mass reconstruction methods? 
There are some numerical indications that the curl-free line- 
integral approaches are close to the KS reconstruction, (Seitz 
& Schneider 1996, Lombardi & Bertin, 1998b) and therefore 
close to the ML as well. For the newly developed regulariza- 
tion methods (Seitz et al. 1998, Bridle et al. 1998) it proba- 
bly not the case. The noise is reduced in these methods, in 
an optimal way that real structures should still emerge from 
the noise. The main advantage of the regularization meth- 
ods is their adaptative spatial resolution depending on the 
local signal-to-noise. However the regularization parameter 
is arbitrary and can potentially kill the small scale fiuctu- 
ations that we are looking for, so the analysis done here 
should not apply. It is therefore valuable to make a study 
of the noise properties for all the reconstruction methods in 
order to search for a better use of the noisy regions of the 
mass maps and to compare their ability to detect the mass 
overdensities. 



7 ACKNOWLEDGMENTS 

I thank R. Scoccimarro, Y. Mellier and T. Erben for a care- 
full reading of the manuscript. 



REFERENCES 

Bardeen, J. M., Bond, J.R., Kaiser, N., Szalay, A.S., 1986, ApJ, 
304, 15 

Bartelmann, M., Narayan, R., Seitz, S., Schneider P., 1996, ApJ, 
464, L115 

Bond, J.R., Efstathiou, G., 1987, MNRAS, 226, 655 (BE) 
Bridle, S.L., Hobson, M.P., Lasenby, A.N., Saunders, R., 1998, 

MNRAS, 299, 895 
Connolly A.J., Szalay. A.S.. Brunn er. R.J., 1998, ApJ, 499, L125 
Erben, T et al., istro-ph/9907134 

Hattori, M., 1997, "Cosmic Chemical Evolution", lAU 187, Ky- 
oto, Japan 



astro-ph/9901191 



Jain, B. Seljak, U., White, S., 
Kaiser, N. 1995, ApJ, 439,L1 
Kaiser, N. 1998, ApJ 498, 26 
Kaiser, N., Squir es. G.. 1993. AdJ 4 04. 441 



Kaiser, N. et al., |astro-ph/9809268 

Le Fevre, O., Hudon, D., Lilly, S.J., Crampton, D., Hammer, F. 

Tresse, L., 1996, ApJ, 461, 534 
Lombardi, M., Bertin, G., 1998a, A& A, 335, 1 
Lombardi, M., Bertin, G., 1998b, A& A, 330, 791 
Lombardi, M., Bertin, G., 1999, A& A, 348, 38 
Oukbir, J., Blanchard, A., 1997, A& A, 317, 1 
Peebles, P.J.E., 1980, The large scale structure of the Universe. 

Princeton 

Reblinsky, K., Bartelmann, M., 1999, A& A, 345, 1 



Reblinskv. K.. Kruse, G., Jain, B., Schneider, P., astro 



ph/9907250 



Schramm, T., Kayser, R., 1995, A& A, 299, 1 
Schneider, P., 1996, MNRAS, 283, 837 



Schneider, P., Seitz, C, 1995, A& A, 294, 411 

Seitz, S., Schneider, P., 1996, A& A, 305, 383 

Seitz, C, Schneider, P., 1997, A& A, 318, 687 

Seitz, S., Schneider, P., Bartclmann, M. 1998, A& A, 337, 325 

Squires, G., Kaiser N., 1996, ApJ 473, 65 

Szapudi, I., Colombi, S., 1996, ApJ, 470, 131 

Van Wacrbekc, L., Bcrnardcau, F., McUier, Y., 1999, A& A, 342, 
15 

Wilson, G., Gole, S., Frenk, C.S., 1996, MNRAS, 280, 199 



